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1.   Introduction 


Sewall  [  4  ] ,  Waugh  [  6  ],  and  Srinivasan  [  5  ]  discuss  least  squares 
estimators  of  a  very  particular  set  of  two  structural  equations  with  no  exoge- 
nous variables  and  the  net  result  of  their  discussion  is  that  such  estimators 
are  unbiased  predictors  of  one  dependent  variable  given  the  other.   Our  pur- 
pose here  is  to  generalize  Sewall 's  main  result  in  several  direction.^,  and 
to  provide  some  ancillary  facts  about  the  structural  equations  defined  in 
(2.1)  below.   Throughout  this  note  we  distinguish  a  random  variable  from  a 
value  assuned  by  it  with  a  tilde;  e.g.  the  random  matrix  ^.   The  symbol  fi 
denotes  the  Kroenecker  direct  product;  e.g.  A  S  B.   And  we  say  that  the 
(m  X  m)  randon  matrix  £  is  "Wishart  with  parameter  set  (H,  n)'' if  4  hes  density 

I  -   ytr  H  £    -jCn-m-l) 

i   ''o^  '='  6.  H  PDS 


I 


n  >  m-1 


otherwise 


where  c   is  a  normalizing  constant, 
o 


2.   Generalization 


Two  questions  immediately  arise:   First,  "^^^hat  is  the  analogue  of  Sewell's 
assertion  a.out  x'^^/x'^x  when  m  >  2 ,  B  is  arbitrary  (m  x  m)  non-singular,  and 
exogenous  variables  are  present?"  And  second,  "UTiat  does  the  answer  to  the 
first  question  imply  .bout  the  conditional  expectation  of  any  single  endoge- 
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nous  variable  given  the  v,  lues  of  all  other  endogenous  variables  and  the  values 
of  the  exogenous  variables?" 

To  make  these  questions  precise,  consider  the  following  system  of  stochastic 
equations : 


iX 


^a)  =  _rz(J)  +  u(J^ 


(2.1) 


where  B  and  r  are  (m  x  m)  and  (m  x  r)  coefficient  matrices,  fixed  for  all  j, 
_z  -^   is  an  (r  X  1)  vector  of  predetermined  variables  and  ^    ^"^  ^    ^^^ 


(m  X  1)  random  vectors.   We  assume  that  {u 


.-.(j) 


j=l,2,...}  is  a  sequence  of  mutually  independent,  identically  Normal  random 


vectors  with  mean  0^  and  PDS  covariance  matrix  £  =  h^  ;  and  B  is  non-singular. 
One  observes  (_^   >  ^^   )  3~^>^  >  •  •  •    but  neither  B,  nor  V   nor  h  is  known  with 


t   t,  t 


certainty.   If  we  partition  a  generic  ^  into  ^  =  (y_    y_^)      with  ^  of  dimensi 


on 


-1 


(p   x    1)  ,    1  _^  p  ^  m-1,    conformably  partition   n^  =  -B        F   into   n  = 
dimension    (p   x   r)  ,    let    f^  =     B        Z   B  and  set 


h 


,    Hj,    of 


.-1 


=11 


%1 


%2 


£22 


r" 


=11  Si2 


=21 


=22 


fill'  Sii  '^^^  ^^  ^  ^^ 


then 


E(iilx2.   ^)    =   h-~   Sll  221^-^2 


h  -^ 


(2.2) 


and 


Var(x^|X2    z) 


ill  =  fOii 


hi   hi  221^ 


(2.3) 
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Now  suppose  we  observe  a  sequence  {(^/^    ,    zT^h  ,    j  =  l,2,...,n}  of  n  _>  r+m 
sample  observations  generated  according  to  (2.2).   Let 

I  =    iL^^K  . .  .  ,l}''h      ,      Z  =  [z(l\...,z(")]   .   z  of  rank  r, 


and 


V=  E  z^J^  z^J)   =  Z  Z' 


Then  it  is  well  known  ([1]  p.  183)  that  given  B,  -T,    and  l,    the  statistics 

P  z'-"    ] 


i  E  i  z'  V"^   and  I    E  ELi^J)  -  P  z^J^]  [^^J^  -  p  -A^^^ 


) 


are  mutually   independent   and   that   P    is   an  unbiased  estimator  of   II   and  e 

J  r  -  =     n-r  = 

an  uni^iased  estimator  of  fil.   Partition  P  and  e   as  follows: 


Si 


22 


P   dim  (p  X  r) 


and 


e  = 


11 


=  21 


Sl2 


;22 


,     e^^  dim  (p  x  p) 


The  analogue  of  the  central  question  of  [  4  ] ,  [  5  ]  ,  and  [  6  ]  may  now 
be  stated:   Is  it  true  that 

Ui[^^'^  \4^^'\^^'''■'h   -   E(P,)  z^"^^^-  E(i^2  122^4"'''^-  ECP^)  U^""''^  ?   (2-^) 
The  answer  is  "yes"  and  this  is  proven  shortly.   We  go  on,  however,  and  do 


considerably  more  in  the  next  section. 


Defining  |i  i  o  ~  In  ~  Sl2  §22  =21  ^^'^  ='[    2   ~   =12  =22'  "^  show  that  the 
joint  likelihood  of  (P,  R,  „»  £-ii  o'  =22^  given  n  and  Q  has  these  properties: 

(1)  P  is  Normal  with  mean  n  and  variance-covariance  matrix 
H  a  V,  and  is  independent  of  R  „,  e   „>  3"*^  £22' 

(2)  e^,  „  is  (p  X  p)  Wishart  with  parameter  set  (H,,,  n-r) 
and  is  independent  of  P,  R,  2>  ^^'^   ^22'  consequently 

(3)  R,  ^  and  l^^   ^^^  jointly  independent  of  P  and  Ij^^^^' 
while  the  conditional  density  of  R,  „  given  z^^   is  Nor- 
mal with  mean  H~  H   and  variance-covariance  matrix 

U   &   e  „  and  the  marginal  density  of  £22  ^^  Wisuart  with 
parameter  set  (2~2»  n+p-r) .   And  £(£22)  =  (n+p-r)fi22- 

(4)  The  density  of  R,  j   unconditional  as  regards  P,  ^22'    ^"*^ 

S1I.2  ^" 

^"ltSi.2-Sn  Si2lS22tii.2-aii  Si2^    -^  Sill  ''-'^ 

where  c"  is  a  normalizing  constant  equal  to 


r  4[n+p-r])  jr      -  2("-^-) 


1 
2 


^(r-p)  Tpr    , 

:2      ^2   r  ihn-r]) 


IL'22'    '=11' 


p^2' 


and 

P 

var  ■ 


r  (a)    =  .1,    r(a4(i-l))-   This  is  the  generalized  multi- 
p      1—1     ^ 
iate  Student  density  first  derived  by  Savage  (  3  ].   Here 


R   „  has  mean  hTJ  H^„  and  variance-covariance  matrix 


K     nas  mectu  ^.-^   ti-^2 

1   0-1  Q  h"''".   (See  Martin  [2  ]) 
n-r-1  =22  ^   "ll" 
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The  above  properties  clearly  imply  that  (2.4)  holds  and  in  addition,  since 

,,   ,-(n+l)|  (n+1)    (n+1.    u-1   .   -,   ,     1 

Var(x^     1^2     »  ^    ^  =  Sn.  imply  that  j^^  z^^^^   is  an  unbiased  estimate 

^  ,.      ,--(n+l)|  (n+1)    (n+1),    ^.  .,   ,  " 

of  •ar(x-,^     1^2    ,  ^     )•   Similrrly,  since  \aT{^^\z)   =  ^^2.  an  unbiased 

estimate  of  this  variance-covariance  matrix  is  — e  Setting  p=l  and  run- 
ning dov.n  the  list  of  assertions  above  answers  the  second  question  posed  at  the 
outset  of  this  section. 

3.   Proofs 


It  is  well  known  (see  [  1  ],  p.  183)  that  the  joinl  likelihood  of  (P,  t_) 
given  y,  and  H,  is 

-  Jtr   H{[P-n]V[P-n]'^}    -  ytr  H  e    ■|(n-r-m-l) 
c  e  •  e  l^l  (3.1) 

where   c    is    a  normalizing   constant.      To   find   the  joint   density   of 

^S»    =1   2'    Sii   2'    =22^    "^  ^^^'^   ^^^   following 

Lemma:      The   Jacobian   J(P,    £   ^  P.    |n    2'    =11   2'    =22''    °^   ^^^   transformation   from 

(P,    e)    to    (P,    R^    2'    ill. 2'    ^22"*    ^^    i=£:22'^" 

Proof:   We  split  the  Jacobian  into  the  product  of  four  transformations  done 

successively:   P  "^  £,  £22  ^  =^2'  =12  "*  5l  2  ^^°   =11  ^  =11.2"  ^^^^^>   ^^S.  ■*  P)  = 

J(e22  ^  422^  "  -^'      ^°"  §1  2  "  =12  =22'  ^^'^  ^   little  algebra  shows  that 

ill  =  ill. 2  ^  Si. 2  i22  Si. 2-   Consequently  J(|^2  "  Si. 2^  =  '^2!^  ^^ 

J(£   ->   e_  )  =  1.   Multiplying  all  four  Jacobians  together  prcves  the  lemma. 

We  can  now  find  the  joint  density  of  (P,  R^_2'  ill. 2'  =22^  ^^°^  ^^^^  °^ 
(P,  |)  by  substituting  in  (3.1)  and  multiplying  by  the  Jacobian  I422 1  •  ^'^^^^ 


rewrite 
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tr  H  e  =  tr  {H^^  e^^  +  H^^  hi   ^  S21  ^12  +  ^22  hl^ 


=  ^^  Sll  S1I.2  ^  ^^  222  hi   ^  ^^  SllfSl.2-SnSl2J|22L^.2-S'l^l2^' 

by  writing  out  tr  H  e,  substituting  using  the  definitions  of  fi„   in  terms  of 

the  H..,  1  <  i,  j  <  2,  and  of  R    and  e,,  „  in  terms  of  tae  e..,  i  <  i,  i  <  2. 
=ij    —      —  =i.z     =il./  =ij    ~  '  -J  — 

The  joint  density  (3.1)  multiplied  by  l^l^  then  gives  that  of  (P,  R   ,  e     , 
I22)  as 

-  ^tr  H{[P-n]V[P-n]'^}     -  |tr  H^.  t,,    „        ^(n-r-m-1) 

I411.2I 

-  |tr  222  £22  I    i(n+P-r-in-l) 
•  e  I422I  ^^-^^ 

-  i^^  SiifSi.2-Sn  Si2]i22fSi.2-Su  S12]'  I    ;h> 

and  the  assertions  1,  2,  and  3  of  the  previous  section  follow  directly. 
To  prove  assertion  4,  rewrite  that  part  of  (3.2)  involving  £^2  ^^ 

-  ¥'  ^  hi    I    i(n+2p-r-m-l) 


wh 


,-1  „    it  „    .„       0-1  „    1   ,  0-1 


to 


-.ere  A  =  [R^^2-  ^11  212^''  SlltSi.2-  Su  Sl2^  ^  =12'    ^^^^8^^-^'^^  "^^^   '^^^P^^' 
P,  £     ,  and  e   ,  and  use  the  fact  that  the  constant  that  normalizes  (3.3)  is 

^(n+p-r) 
c*|a|        ,  with  c*  a  function  of  m,  n,  r,  and  p  aXone. 

Using  a  -.ell  known  determinental  identity  we  may  now  write  the  complete 
density  as 
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-1  -^  t  ^    -  7(n+p-r) 

c"|  [Ri^2-  ill  ^U^hl^h.l-   Sll  S12]  ^  Siil 


where  c"  is  as  defined  in  assertion  4  of  section  2. 
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6^1 


